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ABSTRACT 

The physical parameters of Herbig-Haro jets are usually determined from emission line ratios, ob- 
tained from spectroscopy or narrow band imaging, assuming that the emitting region is homogeneous 
along the line of sight. Under the more general hypothesis of axisymmctry, we apply tomographic 
reconstruction techniques to the analysis of Herbig-Haro jets. We use data of the HH30 jet taken by 
lHartigan fe Morsel (|2007f ) with the Hubble space telescope using the slitless spectroscopy technique. 
Using a non-parametric Tikhonov regularization technique, we determine the volumetric emission line 
intensities of the [S II] AA6716,6731, [O I]A6300 and [N II]A6583 forbidden emission lines. From our 
tomographic analysis of the corresponding line ratios, we produce "three-dimensional" images of the 
physical parameters. The reconstructed density, temperature and ionization fraction present much 
steeper profiles than those inferred using the assumption of homogeneity. Our technique reveals that 
the reconstructed jet is much more collimated than the observed one close to the source (a width 
~ 5 AU vs. ~ 20 AU at a distance of 10 AU from the star), while they have similar widths at larger 
distances. In addition, our results show a much more fragmented and irregular jet structure than 
the classical analysis, suggesting that the the ejection history of the jet from the star-disk system 
has a shorter timescalc component (~ some months) superimposed on a longer, previously observed 
timescale (of a few years). Finally, we discuss the possible application of the same technique to other 
stellar jets and planetary nebulae. 

Subject headings: ISM: Herbig-Haro objects ISM: jets and outflows - Techniques: image processing 
- Methods: data analysis - Stars: pre-main sequence 



1. INTRODUCTION 

Herbig-Haro (HH) jets appear as narrow and well col- 
limated (~ a few degrees) knotty beams of atomic/ionic 
and/or molecular gas moving away from the source 
with supersonic velocit ies (~ 100 — 300 km s _1 ) (e.g. 
IReipurth fc B£dl\l[200l . 

While considerable progress has been obtained in the 
understanding of the behavior of HH jets, long stand- 
ing problems remain unsolved, including the origin of 
the s tellar jets themselves (e.g. lFerreiralll997l : iShu et ahl 
120001 ) . the importance o f jets in solving t he angular mo- 
mentum problem (e.g. iHartmaiml l2009fl , the orig in of 
knots along stellar jets (e.g. iDe Colle et alJl2008bT) . the 
feedback of jets on the core and the star formation ef- 
ficiency, and the interaction of jets with the ambient 
medium eventually leading to the generation of t urbu- 
lence in molecular clouds (e.g. iCarroll et all 120091 ) . To 
properly understand these phenomena, it is crucial to 
accurately determine the physical conditions of the jet 
plasma. 

The proper determination of the electron density n e , 
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temperature T e , and hydrogen ionization fraction ih 
gives strong insights, for instance, on the nature of the 
heating source of the jet material, on the mechanism that 
leads to the production of the knots and on the inter- 
action of the outflowing material with the surrounding 
environment. The determination of the total density «h 
(= n c /xii) is important in order to derive the mass flux, 
and more significantly, to understand the role played by 
jets in the removal of angular momentum from the star 
and the circumstellar disk. Furthermore, different jet 
ejection models predict different dependencies of wh as a 
function of the distance fr om the jet axis (e.g. IShu et all 
fl99llLerv fc Frank) l2l)00h . 

Two approaches are used to infer the outflow ex- 
citation conditions from the observations. The phys- 
ical parameters can be cal culated directly from the 
observed line rat i os (e.g. | Bacciotti fc Eisloffell 119991: 
Nisini et all 120051 iPodio et all 12003 lHartigan fc Morsel 
20071 hereafter HM07). Alternatively, the observations 
can be compared with predictions of emission line inten- 
sities from models of different complexit y, ranging from 
"simple" parallel sho ck models (e. g. | Hartigan et all 
19871 lHartigan et allll994 lLavallev-Fouquet et al.ll2000L 
Pesenti et all 120031) to complex mu lti-species numerical 
simulations (e.g. iRaga et all 120071 ). All these studies 
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agree that the jet is denser close to the source, with tem- 
peratures and ionization fractions which, although peak- 
ing at the positions of the knots, tend to decrease for 
larger distances from the source. 

While previous studies have focused on the variation of 
the physical conditions along the jet axis, recent imaging 
of resolved jets has been performed ( e.g. iBacciotti et all 
20001 [200lTCoffev et al.N200l l2007t IWoitas et alJl200a 
Beck et alll2007l HM07). It has been observed in these 



studies that the electron and total densities peak on the 
axis of the jet and decrease away from the axis (within 
e.g. a factor of 2 in HM07). The temperature and ion- 
ization fraction cross-sections are found to be more ho- 
mogeneous. 

We note, however, that the above-mentioned find- 
ings were obtained assuming that stellar jets are homo- 
geneous along the lin e of sight. In a previous paper 
(|De Colle et al.l l2008al hereafter Paper I) we discussed 
in detail the error in the determination of the physical 
parameters introduced when assuming a jet as homoge- 
neous along the line of sight. In Paper I we proposed 
the application of tomographic techniques to determine 
the volumetric physical parameters, showing an exam- 
ple of the applic ation of a parametric "multi- Gaussian" 
(Bendinelli 1991) method to reconstruct the cross-section 
of the HH30 jet at one particular position along the jet. 
In this paper, we extend the previous results by recon- 
structing the full two dimensional structure of the HH30 
jet by implementing a non-parametric Tikhonov recon- 
struction technique. 

In particular, Section 2 explains in detail the problems 
related to the inversion of the physical parameters, and 
the methods employed in the inversion problem. Section 
3 presents the application of the reconstruction process 
to the HH30 jet. Section 4 and 5 present a discussion of 
the results and our conclusions. 

2. THE DETERMINATION OF THE JET STRUCTURE 

2.1. General considerations 

We consider a jet with known physical parameters 
(hereafter, we refer to the jet electron density n e , elec- 
tron temperature T e and ionization fraction xh as the 
physical parameters of the jet). To determine synthetic 
cmissivitics from the physical parameters it is necess ary 



pa: 

to follow a series of steps (e.g. iDougados et aLl l^OlO) 

1. The emission coefficients are calculated, assuming 
statistical equilibrium and determining the popu- 
lation of each level for each considered atom/ion. 

2. The emission coefficients are integrated along the 
line of sight. 

3. The synthetic emissivity are convolved with the 
seeing (i.e., blurring and scintillation of the jet in 
the Earth's atmosphere), the response of the de- 
tector, the sampling design and the instrumental 
point spread function (PSF). 

The inverse process consists of deriving the emission 
coefficients from the observed emission lines, and deter- 
mining the physical parameter from the ratio of the emis- 
sion coefficients. Actually, the observed emission lines 
always have an associated intrinsic noise. As a conse- 
quence, deconvolution and inversion processes are in gen- 



eral "ill-conditioned" problems (i.e. similar initial con- 
ditions lead to very different deconvolved/inverted so- 
lutions). The level of ill-conditioning depends on the 
particular form of the convolution/integration operator 
("the kernel"), and on the level of noise in the data (e.g. 
iCraig fc BrowrJll986t lBrownlll995h . 

The convolution also produces a loss of information, 
that again depends on the form of the kernel. In fact, 
the convolution process may be seen as a low-pass fil- 
ter. All the structures with size smaller than the char- 
acteristic size of the kernel (the PSF/seeing) will not be 
reconstructed properly when deconvolving the image. 

In addition, the determination of the physical param- 
eters from the emission coefficients depends on several 
assumptions and uncertainties as, for instance, those as- 
sociated wffhtheratiosof the populations of the emitting 
species (jPodio et al l [200fih . 

Finally, to solve the tomography problem, it is nec- 
essary to make some assumptions about the geometry 
of the problem. In some cases, in stellar jets, a rea- 
sonable hypothesis is to consider the jet as cylindrically 
symmetric. This hypothesis describes better the three 
dimensional structure of the jet than the commonly used 
assumption of homogeneity along the line of sight. 

2.2. Reconstruction of the Jet Structure Using 
Tomographic Techniques 

Let us consider an axisymmetric jet moving in the 
plane of the sky. In the following, we use a two- 
dimensional Cartesian reference system centered on the 
star-disk system, and defined on the plane of the sky, 
such that the x- and y-axes are perpendicular and par- 
allel (respectively) to the "main" axis of the jet. Addi- 
tionally, we indicate with r the cylindrical distance from 
the jet axis. 

The aim of the inversion process is to determine the 
emission coefficients i(r,y) as a function of the observed 
intensities I(x, y). The problem is further reduced to one 
dimension assuming that the variation of the physical 
parameters and therefore of the emission lines along the 
y-axis (parallel to the direction of propagation of the jets) 
is much smaller than the variation in the direction across 
the jet (i.e., along the r, x-axis). 

With these approximations, the observed emission line 
intensity cross section I(x) is related to the emission co- 
efficient i(r) by the Abel transform: 



I(x) 
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VT 
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where R is the jet radius, equation ([T]) represents an 
approximation of the relation between observed and vol- 
umetric emission coefficients (see Appendix A). 

As mentioned previously, the inversion of equation ([I} 
constitutes an ill-posed problem. The reason for that is 
easily seen looking at the analytical solution of the Abel 
transform, given by: 
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Decomposing the data into a continuous noise-free and 
an oscillatory noisy part (expanded as a function of the 
frequencies ui n ) I = Iq + S J2 n sinw n x, the derivative 
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is: dl/dx = dlo/dx + S u n sin oj n x, and the high fre- 
quency noise dominates the integral of equation ([2]) . 
Equation ([1]) can be written as 



d = Ax , 



(3) 



where a; is a vector of M elements (a discretization of 
the solution i(r) of the inverse problem), d is a vector 
of N elements representing the data I(xi), and A is a 
matrix with dimension M X N. A is an operator that 
transforms the volumetric into the integrated intensities, 
and is given by a particular discretization of the integral 
present in equation (JTJ) (sec Appendix A). 

The least-squares solution of this linear system of equa- 
tions is obtained minimizing the norm of the residuals: 



\\Ax-d\ 



2 

2 • 



(4) 



where the standard definition of the 2-norm of a vector 
of n elements is given by 



+ x\ + ■ ■ ■ + xl . (5) 

For noisy data, the solution determined minimizing 
the norm of the residuals (eq. [4]), as shown before, is 
dominated by high fre quency noise. A standard solution 
to th is problem (e.g., iCraig fc Brownl IT986I : lAster etahl 
120051 ) is obtained by "regularizing" the data, which in 
the Tikhonov approach consists of solving the following 
damped least squares problem: 



(\\Ax-d\\l+ a 2 \\Lx\\l) 



(6) 



with the condition: 



\Ax- d\\ < 5 



(7) 

where a is the "regularization parameter", L is the 
"Tikhonov matrix", and 5 is the uncertainty (i.e., the 
noise level) in the data. L can be chosen to minimize the 
norm of the solution, or the first or second derivatives of 
the solution (see Appendix A). 

Basically, the Tikhonov method chooses the smoothest 
solution (with respect to the zero-, first-, second-, etc. 
order derivative, depending on the particular form of L 
used) such that the residuals are of the order of the noise 
in the data. The value of a obtained by equation ([7]) 
determines the "optimal" (i.e. the smoothest) solution 
compatible with the data and their associated errors. In 
particular, while a too small value of a produces a so- 
lution dominated by noise, a larger value of a produces 
a too smooth solution. Finally, it is important to note 
that the solution of the inverse problem depends on the 
particular form of the Tikhonov matrix L (that is a free 
parameter of the problem). This is related to the well- 
known fact that, to solve an inverse problem, it is nec- 
essary to introduce a certain amount of information on 
the form of the solution (e.g., on the smoothness of the 
solution or its derivative of order n). 

Once the emission coefficients have been reconstructed, 
the physical parameters are determined from their ratios. 

3. RESULTS 

3.1. The HH30 Data 

The HH30 je t is one of the best s tudied outflows from 
T Tauri stars tiMundt et all 119911: iBurrows etahl 1 1996 : 
iRav et all [1996]; iBacciotti et all 119991 iPetv et ahl I2OO6 ; 
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Fig. 1. — Position of the jet axis 20 relative to the star. The x 
and y-axis are parallel and perpendicular to the main axis of the 
jet, respectively. Negative/ positive x distances correspond to the 
redshifted/blueshifted outflow. 



lAnglada et~aT]|2007l HM07). It is located at a distance 
of 140 pc (jKenvon et al.lll994D . in the Taurus molecular 
cloud. HH30 is nearly on the plane of the sky, within an 
angle of ~ 1°, and it is approximately axisymmetric. Be- 
cause the outflow lies so close to the plane of the sky, the 
flared disk absorbs a large part of the radiation coming 
from the star, allowing the study of the jet at distances 
> 20 AU from the star. The jet is bipolar, but with a 
blueshifted lobe much more luminous than the redshifted 
lobe. 

Furthermore, high resolution, multi-epoch images of 
this jet are available. The data, presented by HM07, 
were obtained by the Space Telescope Imaging Spectro- 
graph (STIS) on HST using a slitless spectroscopy tech- 
nique, during two epochs, in 2000 and 2002. They con- 
sist of two-dimensional images of the jet in a series of 
emission lines, including four bright optical lines: the 
[S II]AA6716,6731 doublet, the [O I]A6300 line and the 
[N II]A6583 line. 

In the images, the cross section of the jet (i.e., in the 
direction perpendicular to the jet axis) is resolved with 
~ 15 — 20 pixels, depending on the distance from the 
central star. As determined by HM07, the PSF FWHM 
is of ~ 4 pixels, with each pixel size corresponding to 
0.025"or 3.55 AU/pixel. As stated by HM07, the S/N is 
very high along the jet knots, with uncertainties below 
5%. 

3.2. Line reconstruction 

The Tikhonov regularization technique has been ap- 
plied to the inversion of the observed data, to reconstruct 
the emission coefficients. 

To invert equation (JT|) it is necessary to know the posi- 
tion xq of the jet axis in the emission line images. As the 
HH30 jet is not perfectly axisymmetric, xq is a slowly 
varying function of y, i.e. xq — xo(y). We determine 
x o(y) by a Gaussian fit to the total line intensity cross- 
section. The resulting Xo is shown in Fig. [TJ 

The blueshifted jet be nds at distances > 40 AU, prob- 
ably due to precession ((Burrows et al.|[l996[) . while the 
redshifted jet (much more noisy due to the lower S/N 
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Fig. 2. — Observed (left) and reconstructed (right) emission line 
images (arbitrary units). The images refers to the sum of all emis- 
sion lines (reconstructed independently). The central source is on 
the top of the Figure. 

of the emission line images) it is not aligned with its 
blueshifted counterpart, forming an angle of ps 2.5°. It 
is unclear whether this asymmetry is due to inhomo- 
geneities of the ambient medium or is a property of the 
ejection mechanism, although the fact that it starts close 
to the source seems to favor the latter explanation. 

An asymmetry of the ejection mechanism itself could 
be related to an asymmetry in the stellar magnetic field. 
However, explor a tory three-dimensional simulations by 
iRomanova et al.l (|2009j ) of a stellar dipole magnetic field 
misaligned with the rotation axis of the star and the disk 
produce winds that are axisymmetric about the rota- 
tional axis. More co mplex magnetic field co nfigurations 
have been studied by lLovelace et al.l (|2010| ) . These au- 
thors used two-dimensional axisymmetric simulations to 
study the case of a rotating star with a large quadrupole 
component, showing that intrinsically asymmetric jets 
(with respect to the equatorial plane) can be generated 
in this case. Therefore, a possible origin for the asymme- 
try could be a large quadrupole component of the stellar 
magnetic field, misaligned with the star + disk rotation 
axis. 

Consistently with the axisymmetry hypothesis, pixels 
symmetric with respect to xq are combined to produce 
symmetric emission-line cross-sections. The emission co- 
efficients are then determined inverting the Abel trans- 
form (eq. [T]) for each emission line. The resulting emis- 
sion line images (Fig. [2]) are sharper than the original, 
while still conserving the main features of the jet (e.g. 
the position and number of knots). 

However, sharper emission lines do not necessarily im- 
ply different line ratios and therefore different physical 




[Nil] 6583/ [SII] 6716+6731 ori 

Fig. 3. — Reconstructed vs. original [N II] /[S II] ratio (upper 
panel) and [S II] 6716/ [S II] 6731 ratio (lower panel). Only pixels 
with emission lines with S/N > 10 arc included. The dashed line 
shows the curve y = x. 



parameters. In fact (see Paper I), in the case of an ob- 
served emission line with a gaussian dependence on x: 
I{x) cx exp(— (x/a) 2 ), the emission coefficients scale as 
i(r) cx l/<7 cxp(— (r/a) 2 ). The relation between the re- 
constructed and original line ratios is therefore 



The [S II] 6716/[S II] 6731 and [N II] 6583/[S II] 
6716+6731 reconstructed ratios are shown in fig. [3] as 
a function of the observed ratios. 

The [S II] 6716/ [S II] 6731 reconstructed ratio is sys- 
tematically lower than the original one, with a larger ef- 
fect for lower values of the ratio. As this ratio is a mono- 
tonic decreasing function of n e , this result implies that 
the inversion process produces in general larger electron 
density, with a larger increase for densities close to the 
high density limit. The large spread in the data is not 
related to errors in the tomographic inversion process. 
The inversion of a pixel with a certain [S II] ratio is not 
unique, as it depends on the emission line cross-sections 
and not only on its own value. 

The [N II] 6583/ [S II] 6716+6731 ratio presents a more 
complex behavior. It depends on all the physical pa- 
rameters (with a stronger dependence on xh), and it is 
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Fig. 4. — Electron density for the first and second epoch for the 
original (left panels) and reconstructed emission line images (right 
panels). The x-axis represents the direction perpendicular to the 
jet axis (with a size of 60 AU), while the y-axis represents the 
direction parallel to the jet axis (500 AU). The source is located at 
the top of the Figure. 

a good shock tracer (HM07), as it becomes higher for 
denser, hotter and more ionized gas. Fig. [3] show that 
this reconstructed ratio increases (with respect to the 
original) for high values, while slightly decreasing for low 
values, indicating that shocks will be more evident in the 
reconstructed jet. 

3.3. The Determination of the Physical Parameters 

To determine n e , T e and xn we use the BE method 
(jBacciotti fc Eisloffell[l999| ). In this method, Nitrogen 
and Oxygen ionization fractions are determined assum- 
ing charge exchange equilibrium with Hydrogen, while 
Sulfur is assumed to be completely ionized. From four 
emission lines, three independent ratios are formed. Be- 
ing the ratios just function of n e , T c and xjj (fixing the 
relative population of S, O, N), these parameters can be 
determined easily by any root-finding algorithm. 

This procedure has been improved recently by HM07, 
including all the available ratios. HM07 showed that the 
best estimation of the physical parameters is the one ob- 
tained by the minimization of the least square problem 



v 

E 

k=l 



(r k -m k y 



(9) 



where = mfe(n c , T c , ih) is the model prediction, r k is 
the natural logarithm of the observed fc— ratio, p is the 
total number of ratios considered, and o~k is the error in 
the fc— observed ratio. 

The electron density, temperature and ionization frac- 
tion obtained from the observed and inverted ratios are 
shown in fig. [4][6] for the two epochs, while the derived 
total hydrogen density is shown in fig. [7] 
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Fig. 5. — The same as in Fig. [4] but for the electron temperature. 
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Fig. 6. — The same as in Fig. [4] but for the ionization fraction. 



In general, the reconstructed n Q , T and ih cross sec- 
tions are much steeper than the observed cross sections of 
the same parameters. The [S II] ratio is in the high den- 
sity limit close to the source, and the electron densities in 
that region are > 2x 10 4 cm -3 . This region can be iden- 
tified as the saturated region in Fig. @] While both the 
observed and the reconstructed electron densities drop 
with the distance from the source, the detailed behav- 
ior is quite different. In fact, the reconstructed electron 
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Fig. 7. — The same as in Fig. [4] but for the hydrogen density 
(n H /n e ). 



density presents a series of stronger on-axis peaks with 
n c > 10 4 cm -3 extending along all the jet. In those re- 
gions, the ionization fraction also increases, while the 
temperature tends to present on-axis valleys. Corre- 
sponding to lower values of T c on the axis, the jet presents 
regions of T c > 10 4 K at m 20 AU from the jet axis. The 
inversion process produces total densities (njj = n e /xn) 
that are also larger on the jet axis with respect to the 
jet edge, although this difference is lower than the corre- 
sponding difference in the electron density. While the n c 
increases by a factor of ~ 2 (as in HM07) from the edge 
to the axis of the jet for the observed data, the recon- 
structed n e increases by a larger factor, of order of ~ 5-10 
depending on the position along the jet. In general, high 
electron density regions are more strongly affected by the 
reconstruction process than low density regions. This is 
direct consequence of the stronger dependence on density 
of the [S II] ratio close to the high-density limit (see also 
Fig. 8 of Paper I). 

The ionization fraction (both observed and recon- 
structed) seems to increase on the edge of the jet up 
to values > 0.5. On the edge of the jet, the S/N becomes 
low and the error in the determination of the physical 
parameters becomes dominant, so this result is probably 
affected by numerical uncertainties. Interestingly, the in- 
crease in the ionization fraction on the edges of the jet is 
a result expected in shock models, as pressure gradients 
in the post-shock region move the ionized gas laterally in 
the cocoon of the jet, where a region of low density, high 
ionization fraction is created. The ionized gas present in 
that region recombines on a timescale much larger than 
on the axis of the jet, where the density is orders of mag- 
nitude larger. 

3.4. FWHM 



The degree of collimation of the jet can be esti- 
mated by the FWHM. Typical T-Tauri star jets, includ- 
ing HH30, have an opening angle of a few degrees at 
~ 100 AU from the so urce, consistent with magnetic 
self-collimation models (jRav et all [2007). HM07 mea- 
sured a deconvolved FWHM of « 10 AU close to the 
HH30 source, with an opening angle of 2.6°. HM07 also 
showed that there is no clear correlation between the jet 
collimation/ width and the presence of a knot. 

We have calculated the FWHM before and after the in- 
version (Fig. HJl for the sum of the [S II], [O I] and [N II] 
emission lines. We have also deconvolved the FWHM de- 
termined from the original and the reconstructed data, 
assuming a FWHM PSF of 4 pixels (HM07). The decon- 
volution of the reconstructed data gives just an approxi- 
mated estimate of the real FWHM, as the deconvolution 
should be applied to the data befo re the reconstruction 
is made and not after (see eq. |Al| ). 

Fig. [5]shows that results similar to HM07 are obtained 
for the FWHM calculated from the observed and the de- 
convolved images. Also, Fig. [8] shows that the inversion 
and the deconvolution produce a comparable effect on 
the FWHM, which has a similar amplitude at small dis- 
tances from the source, going from a width of ~ 15 AU 
at 20 AU to 30 AU at a distance of 500 AU from the 
source. At larger distances the FWHM computed from 
the image inverted by tomography is slightly lower than 
that obtained from the deconvolved image. The value 
of the FWHM determined by deconvolving the recon- 
structed image is much smaller, going for instance from 
5 AU at 30 AU to 25 AU at 500 AU for the first epoch. 
The opening angle is 2.4°, similar to the one determined 
by HM07. 

The reconstruction process also shows some correlation 
between the knots and the width of the jet. In fact, the 
first two knots in the first epoch and the first knot in the 
second epoch (fig. [5]) clearly correspond to a drop in the 
FWHM, while the rest of the knots do not present any 
correlation. That seems to imply that the knot material 
is actually ejected more collimatcd from the star-disk 
system than the material outside the knots. The reason 
for the loss of correlation seems therefore to be related 
to the presence of collisions between knots. 

3.5. Mass Loss 

Assuming a uniform velocity cross-section of 200 km 
s _1 , and using the density stratification coming from the 
inversion process, we have determined the mass-flux by 
integrating dM oc rvndr. The resulting mass-fluxes com- 
puted from the original and the reconstructed densities 
are similar. They rapidly drop from 10 -8 Mq yr _1 close 
to the source to 10 -7 M Q yr" 1 at 200 A U, that are typi- 
cal values for jets from T-Tauri star s (e.g. JHartigan et al.1 
119951 iLavalfev-Foucmet et al.ll2000D . Interestingly, a very 
different result is obtained assuming a power-law pro- 
file for the velocity. If, for instance, v oc r~ v , values 
°f V ~ 1 gives mass fluxes < 10 ~ 8 M Q yr _1 close to 
the source, and approximately one order of magnitude 
lower at at distance of > 200 AU from the source. These 
values of the mass-flux are much less than those ob- 
tained by using other method s fe.g.. lHartigan et aLlll995l 
lLavallev-Fouquet et al.|[2000l) . implying apparently that 
the velocity cross-section has to be nearly flat (i.e. with 
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Fig. 8. — Upper panels: Emission line flux (in arbitrary units) as a function of the distance from the star, for the I (left) and II epoch 
(right). The flux is obtained adding the [S II], [N II] and [O I]. Lower panels: FWHM as a function of distance from the star, for the 
integrated emission line intensities (shown in the upper panels of the figure). The different curves correspond to the observed (dash-dotted 
line), reconstructed (dotted), deconvolved (dashed) and reconstructed and deconvolved (continuous) data. The deconvolution has been 
done assuming a FWHM PSF of 4 pixels. 



a velocity gradient with r\ < 0.5) across the jet cross- 
section. Alternatively, these results can imply that the 
on-axis total density value is much larger than the one 
derived by the reconstruction process. 

3.6. Errors 

There are three main sources of error that affect the 
results of this work: errors due to the hypothesis of ax- 
isymmetery, errors related to the assumptions of the BE 
method, and errors due to the PSF convolution. They 
are discussed in some detail in the following. 

The first source of errors is related to the hypothesis of 
axisymmetry. To understand how asymmetries present 
in the HH30 jet affect our results, we have determined 
the physical parameters for both sides of the jet inde- 
pendently (see £|3.2[) . The physical parameters obtained 
differ from those shown by less than 30%, without any 
systematic effect present. 

The second source of errors is related to the approxi- 
mations used in the BE method. In particular, the main 
errors in the BE method are r elated to the uncer tainty 
in the elemental abundances. IPodio et al.l (|2006[ ) have 
shown that the choice of elemental abundances affects 
the derivation of the physical parameters up to 40% for 
x e and 25% for T e . On the other side, the determination 
of n e , coming from the ratio of sulfur lines, is not affected 
by the choice of the relative abundances. Furthermore, 
the errors in the determination of the electron density 
are larger close to the star, where the [S II] 6716/6731 
ratio is in the high density limit. 

The third source of errors is related to the PSF. In 
this work we have not deconvolved the data with the 
PSF, for two reasons. First, the exact form for the PSF 
is not available. Second, and more important, the mathe- 
matical problem becomes much mo re complex if the PSF 
deconvolution is included (eq. |A1| ) . As we have seen in 
the determination of the FWHM, the effect of the PSF 
is of the same order as the effect produced by the in- 
version process itself close to the source, and is of lower 



importance at larger distance from the source. Anyway, 
the inversion process itself is limited by the presence of 
noise in the data. A large level of noise implies the need 
to use a larger regularization parameter. In other words, 
the presence of a certain level of noise implies a loss of 
information, that cannot be recovered in the inversion 
process. However, the PSF deconvolution would produce 
effects in the same direction as the inversion process, re- 
sulting in steeper profiles, with larger densities and ion- 
ization fraction, and lower temperatures on-axis. The 
results of this paper remain still valid at least qualita- 
tively, and represent a huge improvements with respect 
to previous works in which neither the PSF deconvolu- 
tion nor the line-of-sight integration were considered in 
the determination of the physical parameters. 

4. DISCUSSION 

The main results of this paper are that the inverted 
electron density is in general larger, and that it presents 
a series of sub-structures not clearly evident in the n e 
determined directly from the observed data. The same 
effect is evident in the total density, although in lower 
measure, as the increase in electron density corresponds 
to a (lower) increase in the ionization fraction. 

This result has important implications for our under- 
standing of the origin of HH objects. The presence of 
knots in stellar jets seems to be related to the generation 
of internal working surfaces due to supersoni c variations 
of the ejection velocity at the base of the jet (jRaga et al.l 
I199CA . On the other side, the morphology of the HH30 jet 
is not easily modeled in terms of variable hydrodynam- 
ics jet models, because a periodic velocity ejection from 
the central source is more likely to produce knots with 
"bow-shock-like" mo rphology (e.g., as seen in HH 111, 
IReipurth et al.lll996f). A possible ex planation has been 
proposed by iDe Colle fc Ragal (|200l . who showed that 
axially elongated structures such as the ones observed 
in the HH30 jet can be produced by pinch instability 
driven by strong toroidal magnetic fields (possibly co- 
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existing with the presence of internal working surfaces). 
The results of this paper suggest a possible alternative 
explanation. In fact, the presence of sub-structures in the 
elongated knots indicate that a velocity variation with 
lower timcscalcs, maybe of chaotic nature, can be oper- 
ating, superimposed to the large scale velocity variation 
responsible for the creation of knots. The scenario of 
chaotic ejection of material forming the finally visible jet 
knots h as recently been ex plor ed in numerical sim ula- 
tions bv lYirak et all <|2009Q and lBonito et all (|2010( ). 

The reconstructed jet presents regions of lower temper- 
ature and higher ionization fraction. In particular, lower 
temperature on-axis regions are qualitatively consistent 
with shock models, due to the stronger cooling that one 
would expect in higher density regions. 

The same technique could in principle be applied to 
other existing data, although other jets have in general 
inclination angles different from zero and, as a conse- 
quence, the int erpretation of the results would be more 
comp lex fe. g. iBacciotti et al.ll2000l 120021: [Coffey et al " 



2004 : IWoitas et all 120051 ICoffev et all 120071: IBeck et al 
20071 HM07). 



The HH30 data presented in HM07 have a very high 
S/N level. While the same technique can in principle be 
applied to data with higher uncertainties, our ability to 
understand the structure of the jet is limited by the noise 
level of the observations, as the regularization techniques 
give smoother solutions for data with lower S/N. 

Observations obtained with integral field spectro- 
scopes (i.e. three-dimensional spectral data with two- 
dimensional spatial coverage), in particular, can be used 
to extract information also on the dependence of the ve- 
locity on r. In this case, in fact, the relation between the 
observed intensity and the emission coefficient is given 

/(^H2f%^S, (10) 



with 



y/ r 2 



j( rjW ) = _*W- e -(*-«(r))V^ 



(11) 



where a v is the velocity dispersion. Given the observed 
line intensity I(x,v), i(r,v) can be determined from 
equation (|10[) , i(r) can be determined by inverting the 
Abel equation (related to eq. [11] by I(x) = J I(x,v)dv) 
and v(r) can be finally inferred from equation (|11[) . In 
practice, a v is often dominated by the instrumental PSF 
(e.g. 10 — 20 km s" 1 vs. a thermal width of ~ a few 
km s _1 for ions as S II or O I). On the other side, often 
v 7S> <J V (~ 100 km s _1 both across and along the jet) and 
important information about the dynamics of the jet can 
be determined from v(r), as for instance the mass-flux. 
In cases when information on the width of the jet in 



two different lines (i.e. the [S II] doublet) is available, 
equation flSJ can be used (see also eq. [18] of Paper I) 
to have an order of magnitude estimation of the effect of 
projection effects. 

5. CONCLUSIONS 

In this paper we have for the first time applied the to- 
mographic inversion technique to determine, under the 
hypothesis of axisymmetry, the three dimensional struc- 
ture of the HH30 jet. This hypothesis describes better 
the three dimensional structure of the jet than the com- 
monly used assumption of homogeneity along the line of 
sight. 

The main results are summarized as follows: the recon- 
structed electron and total densities show a fragmented 
structure, probably consistent with small timescale (~ 
months) ejection velocity variations. Corresponding to 
the peaks in electron and hydrogen density, the ioniza- 
tion fraction also increases, while the electron tempera- 
ture presents on-axis valleys. The width of the recon- 
structed jet is lower than the corresponding width in- 
ferred directly from the observations. In particular, pro- 
jection and instrumental effects have a similar effect on 
the FWHM. Combining the two effects, the FWHM be- 
comes < 10 AU close to the source. Although our de- 
termination of the mass-flux is limited by the lack of 
information on the velocity radial stratification, we have 
shown that the effect of the density stratification on the 
mass-flux evaluation is negligible. 

The application of the BE technique to the HH30 jet is 
limited to distances from the source > 40 AU, where the 
[S II] 6716/6731 ratio is not in the high density limit, and 
is possible to properly determine the electron density. 

Future observations dedicated to the study of the jet 
close to the source, using line ratios with larger criti- 
cal densities, may allow a direct comparison between the 
results of the inversion and predictions of jet ejection 
models. In this way, it will be possible to compare for 
example the density profiles obtained from the inversion 
with theoretical predictions. For instance, the disk-wind 
and X-wind predict different power-law dependencies of 
the density as a function of the distance from the jet axis, 
while the ma gnetic tow er mo del predicts a sharp outer 
jet edge (e.g. IShu et all 11993 ILerv fc Frankll27)00h . 

We have shown that tomographic reconstruction can 
be used to reconstruct the cross section of the flow, ob- 
taining a description of the three dimensional structure 
of the jet. The use of high resolution observations will 
make it possible to apply this technique to a larger sam- 
ple of objects. 



We thank Pat Hartigan for sharing with us his HST- 
HH30 data. 



APPENDIX 

IMPLEMENTATION OF THE TIKHONOV REGULARIZATION 

Here we present the de tails of our implem e ntation. For a more g eneral description of the Tikhonov regularization 
technique sec for instance lAster et all (|2005f ): ICraig fe Brownl (|1986[ ). 
The observed intensity I{xi,yj) is related to the emission coefficients i(r,y) by the following expression: 



i+1/2 rVj+i/2 
dx I dy 

-i/2 Jyj-1/2 



OO /> + OG 



ip(x — x' , y — y')dx'dy' 



i(r, y)rdr 



x'-x (y) y/r 2 - (x' - X (y)y 



(Al) 
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The first two integrals of equation (|A1[) represent the integration on the pixel centered on (xi,yj) with dimension 
(^1-1/2)^1+1/2) x {yj-i/2i Vj+x/2)- The third and fourth integrals are the convolution with the seeing and the instru- 
mental response ip, and the last term represents the integration along the line of sight, being xo(y) the position of the 
center of the jet (in general a function of the y coordinate). 

To simplify the problem, we neglect the effect of the instrumental response, assuming ip(x—x', y—y') ~ S(x—x', y—y'), 
being S(x) the delta function. Furthermore, we reduce the problem to one dimension assuming i(r,y) « i(r,yj), i.e. 
that the variations in the emission coeffi cien ts along y are much smaller than those along r. 

With these approximations, equation (|A1[) reduces to 



dx 



-1/2 



R 



i(r, yj)rdr 



(A2) 



1 (the exact value of Ayj is not important, as it cancels out 



where Syj = J/j+i/2 — Vj—i/2- We further assume Ayj 
in the line ratios), and we indicate I(xi, yj) = U and i(r, yj) = i(r). 
As shown in £ 12.21 equation (|A2[) is inverted using the Tikhonov regularization technique, i.e 



damped least square problem 



with the condition 



min (|| At - d\\\ + a 2 
d\\ < S 



Lx 



\Ax 



In matrix form, equation (|A3|) can be written as 



(A T A + a 2 L)x = A T d 



solving the following 
(A3) 

(A4) 
(A5) 



where A T is the transpose matrix of A. The solution is obtained inverting this linear system of equations. 

As discussed in £|2.2[ the matrix L can be chosen to minimize the norm of the solution, or its first or second derivative. 
In the first case (zero-order) solutions with low magnitude of x are favored, while in the cases of first- and second- order 
regularization, the magnitude of the gradient and Laplacian are penalized, respectively. Being the solution unknown, 
the choice of the particular form of L is somehow arbitrarily. In this work, we decide to minimize the first derivative 
of the solution as, in our implementation, it produces a less noisy solutions with respect to zero- and second-order 
minimization, and it works better in our tests. The matrix L is therefore given by: 



L 



-1 1 







-1 



-1 1 







1 



To complete the description of the algorithm, we need to specify the form of the matrix A. Let us assume that we 
have N values of If., defined in th e po sitions Xf. (with < k < N — 1, where xq = Axq/2). The external radius is 
defined as 2^+1/2 = R- Equation (|A2[) can be therefore written as: 



r x k+i/2 
I k = 2 J dx 

'1/2 



x k+l/2 



ik{r)rdr 
y/r 1 — x 2 



N 

E 

j=k+i 



ij(r)rdr 
\Jr 2 — x 2 



Taking ik{r) = io,k constant inside each shell, this equation leads to: 



h = io,k ( S 



fe+1/2 
fe+1/2 



S 



fc+1/2 
fc-1/2 



N 

2 io ^{ 



J + l/2 
fe+1/2 



a k+l/2 



J+l/2 
fe-1/2 



cU-1/2 
°fc-l/2 



J=fc+1 



(A6) 



(A7) 



where S, 



•^m \j x2 \ 

x 2 n arcsin x/x n 



x„ arcsin x r 



2/^ 



•c 2 dx 



/x n , and the following relation has been used: 

x* + x" n arcsin x/ x n = S n . The matrix A can be easily deduced from equation (|A7j) . 
We have run a series of tests to verify the accuracy of the algorithms developed in this paper. Fig. |9] shows one of 
these tests. Assuming a certain radial dependence for electron density, temperature and ionization fraction, we have 
computed synthetic emission line intensities. We have then integrated the volumetric emission line along the line of 
sight using equation (|A2|) . and then added a Poisson noise to the data. Finally, we have reconstructed the electron 
density, temperature and ionization fraction using the tomographic technique and the BE method f ^3.3|) . As shown in 
Fig. [9l the resulting physical parameters agree with the original one within 10%. For comparison, we show in Fig. [9] 
also the physical parameter cross-sections obtained assuming the medium as homogeneous along the line of sight (i.e., 
those obtained from the "observed" emission lines). The tests show that also with a 10% noise level in the data our 
implementation of the tomographic inversion gives acceptable results. 
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